Optimum spin-squeezing in Bose-Einstein condensates with particle losses 
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The problem of spin squeezing with a bimodal condensate in presence of particle losses is solved 
analytically by the Monte Carlo wavefunction method. We find the largest obtainable spin squeezing 
as a function of the one-body loss rate, the two-body and three-body rate constants, and the s-wave 
scattering length. 
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Spin squeezed states, first introduced in [ij, general- 
ize to spin operators the idea of squeezing developed in 
quantum optics. In atomic systems effective spins are 
collective variables that can be defined in terms of two 
different internal states of the atoms 2] or two orthogo- 
nal bosonic modes [3] . States with a large coherence be- 
tween the two modes, that is with a large mean value of 
the spin component in the equatorial plane of the Bloch 
sphere, can still differ by their spin fluctuations. For an 
uncorrelated ensemble of atoms, the quantum noise is 
evenly distributed among the spin components orthog- 
onal to the mean spin. However quantum correlations 
can redistribute this noise and reduce the variance of one 
spin quadrature with respect to the uncorrelated case, 
achieving spin squeezing. Besides applications in quan- 
tum communication and quantum information 4], these 
multi-particle entangled states have practical interest in 
atom interferometry, and high precision spectroscopy [H| 
where they could be used to beat the standard quantum 
limit already reached in atomic clocks . 

Different techniques to create spin squeezed states in 
atomic systems have been proposed and successfully real- 
ized experimentally including transfer of squeezing from 
light to matter 7] and quantum non demolition mea- 
surements of the atomic state Q. To go further, it was 
shown that coherent interactions between cold atoms in 
a bimodal Bose-Einstein condensates Q , can in principle 
provide a huge amount of entanglement and spin squeez- 
ing. It is thus important to determine the ultimate lim- 
itations imposed by decoherence to the maximum spin 
squeezing that can be obtained by this method. Several 
forms of decoherence may be present in the experiment. 
The case of a dephasing perturbation was studied in 0. 
In this work we deal with particle losses, an unavoidable 
source of decoherence in cold atom systems, due e.g. to 
collisions of condensed atoms with the hot background 
gas, and to three-body collisions leading to molecules. 

As shown in [3] , bimodal Bose-Einstein condensates re- 
alize the one-axis twisting model proposed in [l[ to create 
spin squeezing. This exactly solvable model predicts a 
perfect squeezing in the limit of a very large system: for- 
mally — + for ^ oo, A^ being the number of parti- 
cles in the system and the squeezing parameter defined 



in Eq.®. We expect losses to degrade the squeezing Q 
that is < ^l-^i^ for any value of A^. However, 

as [Qgg ^ as A^ — > cx), this inequality does not tell 
us what will be the maximum squeezing in presence of 
losses. In particular the limit limAr^+oo Cwith loss could 
be zero (perfect squeezing), a very small constant, or a 
constant close to one (one meaning no squeezing). We 
show that the second possibility is the correct one. The 
best achievable squeezing is obtained for N +oo, and 
we derive its explicit expression, as a function of the scat- 
tering length and the loss constants Ki, K2, K3. 

We consider two spatially separated symmetric con- 
densates a and b prepared in an initial state with A^ par- 
ticles and a well defined relative phase 
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We assume that (f) = initially. Correspondingly, the x 
component of the collective spin Sx = (a^6 + b'^a)/2 has 
a mean value (Sx) — N/2. Here we assume that no exci- 
tation is created during the preparation process and we 
neglect all the other modes than the condensate modes 
a and b. When expanded over Pock states \Na,Ni,), the 
state ([Ij shows binomial coefficients which, for large A^, 
are peaked around the average number of particles in a 
and b, Na — Nb = N/2. We use this fact to approximate 
the Hamiltonian with its quadratic expansion around Na 
andTVfc [11]: = J^.^^^^ EiN,)+^^,{N,-N,)^^l'^{N,- 
A^e)^ where /i^ is the chemical potential for the e conden- 
sate and fi[ = {dN^fJ-e)]^ - In the symmetric case, we can 
write 



Ho = fia^a 
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where x = /^a/^- The first term in Hq is some function / 
of the total atom number: It commutes with the density 
operator p of the system and can be omitted. 

In presence of one, two and three-body losses, the evo- 
lution of the density operator, in the interaction picture 
with respect to Hq, is ruled by the master equation 
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where p = e^Hot/hp^-iHot/h^ 

similarly for b, 7^™) = ^ / d^r|(/)(r)p™, where if™ is 
the m-body rate constant and 0(r) is the condensate 
wavefunction in one of the two modes. In the Monte 
Carlo wavefunctions approach we define an effective 
Hamiltonian iJgff and the jump operators ji™'' 
3 
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We assume that a small fraction of particles will be lost 
during the evolution so that we can consider x and 7^™) 
(m = 2, 3) as constant parameters of the model. The 
state evolution in a single quantum trajectory is a se- 
quence of random quantum jumps at times tj and non- 
unitary Hamiltonian evolutions of duration Tj : 

■ ■ ■ 4'"^nii)e"*^^""'^'*IVX0)) . (6) 

The expectation value of any observable O is obtained 
by averaging over all possible stochastic realizations, that 
is all kinds, times and number of quantum jumps, each 
trajectory being weighted by its probability 

{d)^j2 [ ^^i'^^^ • • • dh J2 immit)) ■ (7) 

k •^0<*i<*2<-*'.<*{e,,m,} 

We want to calculate spin squeezing. In the considered 
symmetric case with zero initial relative phase, the mean 
spin remains aligned to the x axis (Sx) = {b^o), and the 
spin squeezing is quantified by the parameter [3, Q 
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where Sg = {cose)Sy + {sme)S,, Sy = {a^ - b'' a) / {2i) , 
= {a^a - b%)/2 and TV = a^a b%. The non corre- 
lated limit yields = 1, while < 1 is the mark of an 
entangled state 0,01 • In all our analytic treatments, it 
turns out that AS'f = (TV)/4. This allows to express 
in a simple way: 
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With one-body losses only, the problem is exactly solv- 
able. Following a similar procedure as in [llj . we get 



7^ +X[7sin(xt) -hxcos(xi)]e-''*"'^^"^ 
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72 + 2x[7sin(2xt) + 2x cos(2xt)]e-^* 
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The key points are that (i) H^s is proportional to N so 
it does not affect the state, and (ii) a phase state \4>) is 
changed into a phase state with one particle less after a 
quantum jump, Ca.b\4') ^ \4'1- Xt/'^) where t is the time 
of the jump, the relative phase between the two modes 
simply picking up a random shift =Fx''/2 which reduces 
the squeezing. 

When two and three-body losses are taken into ac- 
count, an analytical result can still be obtained by using 
a constant loss rate approximation 11 1 
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We verified by simulation (see FiglT]) that this is valid for 
the regime we consider, where a small fraction of particles 
is lost at the time at which the best squeezing is achieved. 
In this approximation, the mean number of particles at 
time t is 

(iV) = iV[l-^r(™)<] ; r(™) = (iV/2)'"-im7(") (14) 



where r^"^h is the fraction of lost particles due to m-body 
losses. Spin squeezing is calculated from ^ with 
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where the operator N = {N — da) acts on the functions 



Fp{a) = exp 



sin(m/3xi) 
m(3xt cos™ {/3xt) 



(18) 



and all expressions should be evaluated in a = In Na ■ 

We want to find simple results for the best squeezing 
and the best squeezing time in the large N limit. In 
the absence of losses [l| the best squeezing and the best 
squeezing time in units of 1/x scale as iV~2/3^ then 
set N = and rescale the time as xt = t^'^- We 
expand the results (fT2|) and pStilTp for e <C 1 keeping 
r'"™Vx constant, and we obtain in both cases 
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For equal loss rates r*^™\ the larger m, the more the 
squeezing is affected. Introducing the squeezing ^o(i) in 
the no-loss case, the above result can be written as 
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This shows that (i) the fact that only a small fraction of 
atoms is lost at the best squeezing time does not imply 
that the correction on the squeezing due to losses is small; 
(ii) the more squeezed the state is, the more sensitive to 
the losses. In presence of losses, the best squeezing time 
and the corresponding squeezing are 
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In order to find optimal conditions to produce spin 
squeezing in presence of losses and set the ultimate lim- 
its of this technique, from now on, we assume that the 
number of particles is large enough for the condensates 
to be in the Thomas-Fermi regime so that 
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where uq — ^Jh/Moj is the harmonic oscillator length, 
M is the mass of a particle and Co is the geometric mean 
of the trap frequencies, 
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We first analyze the dependence of squeezing on the ini- 
tial number of particles, separating for clarity one, two 
and three-body losses. FiglT] shows the best squeezing 
C^(ibest) as a function of N when only one kind of losses 
is present. The curve without losses is also shown for 
comparison. According to FigdJ one-body losses do not 
change qualitatively the picture without losses and we 
have 'C^(ibost) oc N~^/'^^ for N ^ 00. In the same limit. 
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FIG. 1: Spin squeezing obtained by a minimization of ^ 
over time, as a function of the initial number of particles, 
without loss of particles (solid line), with one-body losses 
(dashed line), with two-body losses (dotted line), with three- 
body losses (dash-dotted line) respectively. Parameters: a = 
5.32nm, lj = 2% x 200Hz, Ki = O.ls'^ K2 = 2 x IQ-^^mVs 
13], K3 = 18 X IQ-^^mVs. The symbols: crosses (plus) are 
results of a full numerical simulation with 400 Monte Carlo 
realizations for two-body (three-body) losses. 



with two-body losses, 'f^(tbest) is independent of A^. With 
three-body losses, S,^{thest) oc N'^^^^ for A — > cx), imply- 
ing that, for a fixed Co, there is a finite optimum number 
of particles for squeezing. 

We now turn to a full optimization of squeezing over 
Co and A^ in the simultaneous presence of one, two and 
three-body losses. To this end, we note that the square 
brackets in Eq. (l23p is an increasing function of C, we 
can then optimize ^^(ti^ggj) by minimizing C with respect 
to Co. Under the conditions Ki ^ and K3 ^ 0, the 
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It turns out that Cmin is proportional to N and 
^^(ibcst, <^°P*) is a decreasing function of A^. The lower 
bound for reached for A^ = 00 is then 
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In practice, one can choose A^ = A^,, in order to have 
^2 = (1 -f 7y)min^2 g_ jy ~ 10%), and then calculate 
the corresponding optimized frequency aj°P* with (j30|) . 
For a suitable choice of the internal state, in an optical 
trap, the two-body losses can be neglected K2 = 0. One 
can get in this case very simple formulas for the optimized 
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parameters and squeezing. For 77 = 10% 14j : 
17.833 ha 
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We now ask whether we can use a Feshbach resonance 
to change the scattering length (but also K3) to improve 
the squeezing. In Figl2] we plot the squeezing parame- 
ter vs the scattering length a. Predicted values of K3, 
as a function of a, are taken from [l^ for ^^Rb in the 
state \F = l,mF = 1) and Ki = O.Ols"^. We calculate 
^opt -j-jjg number of particles needed for 77 = 10% for 
each point in the curve. The dip giving large squeezing 
corresponds to a strong decrease in K3 around 1003. 5G 
(Jfa ~ 3 X 10"^^m*^/s). Close to the Feshbach resonance 
the squeezing gets worse as K3 increases (even if in the 
figure we do not enter the regime K3 ~ ha'^/m). 

Finally we consider the problem of the survival time of 
a spin squeezed state in presence of one-body losses. We 
imagine that the system evolves in two periods: for t < Ti 
the system is squeezed in presence of interactions (x 7^ 0) , 
one and three-body losses; and for t > Ti the interaction 
is stopped (x = 0), e.g. by opening the trap, and the 
system only experiences one-body losses. As t can be 
arbitrarily long, we use the exact solution for t > Ti 
while for the i < Ti ~ ibost, we use the approximation 
(fn|. Then for t = Ti+T2> Ti: 
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This result shows that the spin squeezing can be kept 
some time after the interactions have been stopped. To 
give an example, for ^'^Rb atoms with bare scattering 
length a = 5.32nm, Ki = 0.01s-\ K2 = 0, K3 = 6 x 
10-''2m*^/s [il], in optimized conditions (I32l)-(IS1) N 
2.8 X 10^ and (I>°p* = 27r x 20.06Hz, = 5.7 x 10"^ is 
reached at Ti — thest = 4.4 x 10~^s, and a large amoimt 
of squeezing ~ 0.01 is still available after Is. 

In conclusion, we found the maximum spin-squeezing 
reachable with cold atoms having a S'^ Hamiltonian, in 
presence of decoherence (losses) unavoidably accompany- 
ing the elastic interaction among atoms. The best squeez- 
ing is reached for an atom number N ~^ 00 and not for a 
finite value of N. This is important for applications such 
as spectroscopy where, apart from the gain due to quan- 
tum correlations among particles (squeezing), one always 
gains in increasing N . 

LKB is a unit of ENS and UPMC associated to CNRS. 
We acknowledge discussions with the atom chip team 




FIG. 2: Spin squeezing ^^(fbGst) optimized with respect to Co 
as a function of the scattering length a, when the magnetic 
field is varied on the left side of the Bo = 1007. 4G Feshbach 
resonance of *^Rb. The inset shows the number of particles 
for each point, calculated for 77 — 10%. We took a{B) = 
abg[l - AB/{B - Bo)] with Sbg = 5.32nm, AB = 0.21G. The 
three-body rate constant Ki{B) is taken from p^. 
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